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We  evaluate  in  this  work  the  effect  of  Gd  substitution  on  the  dc  activation  energy,  Edc,  and  ionic  conduc¬ 
tivity  of  the  pyrochlore-type  gadolinium  zirconate,  P-Gd2Zr207.  Several  compositions  with  the  general 
formulae  Gd2-yLnyZr207  (Ln  =  Er3+,  Y3+,  Dy3+,  Sm3+,  Nd3+  and  La3+)  were  prepared  by  mechanically  milling 
and  firing  stoichiometric  mixtures  of  the  corresponding  elemental  oxides  and  their  electrical  properties 
analyzed  as  a  function  of  frequency  and  temperature  by  using  impedance  spectroscopy.  Whereas  Gd  sub¬ 
stitution  in  P-Gd2Zr207  by  smaller  lanthanides  induces  a  pyrochlore  to  fluorite  phase  transition,  using 
larger  dopant  cations  yields  partially  disordered  pyrochlores.  Despite  of  higher  structural  disorder,  ionic 
conductivity  values  measured  for  the  fluorite-type  materials  are  lower  than  those  observed  for  pyrochlores 
whereas  activation  energies  for  oxygen  migration  in  the  series  decrease  monotonically  as  the  average  size 
of  the  A  cation  increases. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  highly  efficient  and  environ¬ 
mentally  friendly  energy  conversion  devices  with  a  large  variety 
of  potential  applications.  Current  SOFCs  technology  is  based  on 
the  LaxSri_xMn03  (LSM)  cathode,  the  yttria  stabilized  zirconia 
(YSZ)  electrolyte  and  Ni-YSZ  cermet  anode.  As  cells  are  operated 
at  high  temperatures  under  either  oxidising  or  reducing  atmo¬ 
spheres  over  a  long  period  of  time,  different  materials-related 
problems  associated  with  such  aggressive  environments,  limit  their 
long-term  stability.  Thus,  solid  state  reactions  taking  place  at  the 
cathode/electrolyte  interface  cause  the  formation  of  insulating 
phases  such  as  La2Zr207  and  SrZr03,  both  having  a  detrimental 
effect  on  cell  performance.  A  possibility  to  overcome  this  problem 
would  be  using  the  LSM  cathode  but  with  a  different  and  com¬ 
patible  electrolyte.  In  this  context,  oxides  A2B20(1)60(2)  with  the 
pyrochlore-type  crystal  structure,  represent  an  interesting  group  of 
solid  oxide  ion  conducting  materials  because  of  their  high  composi¬ 
tional  diversity,  structural  flexibility  and  intrinsic  concentration  of 
oxygen  vacancies  [1-3].  The  ideal  cubic  pyrochlore  structure  might 
be  derived  from  that  of  an  anion  deficient  fluorite  by  doubling  the 
unit  cell,  removing  one  out  of  every  eight  anions  and  placing  cations 
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and  anions  in  four  fully  occupied  and  non-equivalent  sites:  16c  and 
16d  for  cations  (Wyckoff  notation)  and  48/ and  8 a  for  anions  (see 
[4]  for  a  detailed  description  of  the  pyrochlore  crystal  structure). 
There  is  a  third  position  in  the  unit  cell  available  for  anions,  8 b, 
which  is  systematically  empty  in  ideal  ( ordered )  and  thus,  non- 
conductive  pyrochlores  but  partially  occupied  in  defect  ( disordered ) 
anion  conducting  materials.  According  to  several  computational 
calculations,  the  most  stable  intrinsic  anion  defect  in  pyrochlores 
is  an  oxygen  Frenkel  pair  consisting  of  a  vacant  48 f  position  and 
an  interstitial  oxygen  ion  located  at  the  8 b  site  [5-8].  The  energy 
needed  to  create  this  defect  might  be  significantly  lowered  by  dis¬ 
ordering  the  cation  substructure,  that  is,  by  inducing  homogeneous 
cation  environments  of  the  oxygen  atoms  [1  ].  Thus,  partially  disor¬ 
dered  phases  are  better  ionic  conductors  than  the  fully  disordered 
fluorite  materials  of  the  same  composition  because  of  smaller  acti¬ 
vation  energies  for  migration;  e.g.,  the  Edc  is  significantly  lower  in 
the  defect  pyrochlore-type  gadolinium  zirconate  (P-Gd2Zr207)  than 
in  its  fully  disordered  analogue,  the  fluorite-type  F-Gd2Zr207  [9]. 
Interestingly,  structural  disorder  in  this  type  of  compounds  is  to 
a  great  extent  linked  to  the  A  to  B  size  mismatch  (the  Ra/Kb  size 
ratio)  and  thus,  materials  showing  different  degrees  of  disorder  are 
possible  by  using  the  appropriate  cation  combinations.  Therefore, 
pyrochlores  are  an  excellent  model  system  to  analyze  the  influence 
of  structural  disorder  in  ion  mobility. 

Although  there  is  a  growing  interest  in  this  type  of  materials, 
the  understanding  of  the  effect  of  A  and/or  B  substitutions  on  oxide 
ion  migration  is  inadequate  at  this  point.  Atomic  scale  comput¬ 
ing  simulations  have  predicted  the  activation  energy  for  oxygen 
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migration  in  A3+/B4+  pyrochlores  to  be  strongly  dependent  on  the 
average  B-site  cation  size  and  give  a  secondary  role  to  the  A  cation 
[5].  However  experimental  data  do  not  confirm  clearly  these  pre¬ 
dictions.  Particularly,  the  effect  of  A-site  homovalent  substitutions 
on  the  activation  energy  and  conductivity  of  pyrochlores,  has  not 
been  completely  understood  and  studies  have  mostly  focused  in  the 
A2(Bi_yBV)207  systems,  with  two  tetravalent  cations  (B,  B'  =  Sn4+, 
Ti4+,  Zr4+)  simultaneously  occupying  the  hexa-coordinated  16d 
site  and  a  single  trivalent  cation  (A  =  Y3+  or  Ln3+)  on  the  octa- 
coordinated  16c  site  [1].  The  aim  of  the  present  contribution  is  to 
study  the  role  played  by  the  average  A  cation  size  on  the  electri¬ 
cal  properties  of  P-Gd2Zr207,  a  well-known  high-temperature  solid 
oxide  ion  conductor.  As  the  size  of  the  B-cation  would  be  constant 
when  replacing  Gd  by  another  homovalent  lanthanide,  electrical 
properties  would  be  mostly  affected  by  the  average  A  cation  size 
and  its  influence  on  structural  ordering/disordering.  As  we  have 
already  shown  [10,11],  La  incorporation  into  Gd2Zr207  decreases 
the  dc  activation  energy  and  even  increases  its  ionic  conductiv¬ 
ity  at  low  temperatures.  We  will  now  extend  this  study  to  some 
other  trivalent  lanthanides  (Er,  Y,  Dy,  Nd  and  Sm)  and  compare  these 
results  with  those  previously  obtained  with  La. 

2.  Experimental 

A  set  of  Gd2_yLnyZr207  materials  with  Ln  =  Er3+,  Y3+,  Dy3+,  Sm3+, 
Nd3+  and  La3+  (see  Table  1  for  a  complete  list)  were  prepared  by 
mechanically  milling  stoichiometric  mixtures  of  high  purity  ele¬ 
mental  oxides  Gd203,  Zr02  and  Ln203  (Ln  =  La,  Nd,  Sm,  Dy,  Y,  Er), 
as  described  in  previous  works  carried  out  in  this  group  (see  for 
example  [12,13]).  Powders  were  mixed  and  dry  milled  together  in 
a  planetary  ball  mill  by  using  stabilized  zirconia  vials  and  balls 
and  a  moderate  rotating  disc  speed  (350  rpm).  Phase  evolution 
of  milled  samples  was  followed  by  using  X-ray  power  diffraction 
(XRD)  in  a  Philips  X’Pert  diffractometer  using  Ni-filtered  Cu  Ka  radi¬ 
ation  (A  =  1.5418  A).  The  mechanically  induced  chemical  reactions 
were  considered  completed  when  no  traces  of  the  starting  reagents 
were  observed  by  this  technique.  Impedance  spectroscopy  mea¬ 
surements  were  carried  out  on  pellets  (10  mm  diameter  and  ~2  mm 
thickness)  prepared  by  uniaxial  pressing  of  the  fine  powders 
obtained  by  milling  (sintering  temperature  =  1500  °C).  Colloidal  Pt- 
paint  was  coated  on  both  sides  of  the  pellets  to  serve  as  blocking 
electrodes.  Impedance  data  were  recorded  over  the  100  Hz-1  MHz 
frequency  range  as  a  function  of  temperature,  by  using  a  Solartron 
1260  Frequency  Response  Analyzer.  Sintered  pellets  were  also 
examined  by  scanning  electron  microscopy  in  a  Philips  XL30  ESEM 
microscope  equipped  with  an  EDAX  Inc.  energy-dispersive  X-ray 
detector  for  microanalysis. 

3.  Results  and  discussion 

Depending  on  the  Ln3+  cation  size,  rare-earth  zirconates 
Ln2Zr207  are  known  to  exist  as  either  pyrochlores  (Ln  =  La  to  Gd) 
or  anion  deficient  fluorites  (Ln  =  Tb-Lu  and  Y3+).  In  addition,  some 

Table  1 

Compositions  of  all  Gd2-yLnyZr207  samples  analyzed  in  this  work.  Ionic  radii 
shown  are  taken  from  [16],  assuming  a  common  coordination  number  equal  to  8 
(Red  =  1.053  A). 


Ln3+ 

R(A) 

Chemical  formula 

Er 

1.004 

Gd2_yEryZr207  (y  =  0.2,  0.4,  0.8) 

Y 

1.019 

Gd2-yYyZr207  (y  =  0.2,  0.4) 

Dy 

1.027 

Gd2_yDyyZr207  (y  =  0.2,  0.4,  0.8) 

Sm 

1.079 

Gd2-ySmyZr207  (y  =  0.2,  0.4,  0.8) 

Nd 

1.109 

Gd2-yNdyZr207  (y  =  0.3,  0.6,  0.8) 

La 

1.160 

Gd2_yLayZr207  (y  =  0.2,  0.3,  0.4) 
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Fig.  1.  XRD  patterns  obtained  for  selected  Sm-  or  Dy-containing  Gd2-yLnyZr207 
samples  after  firing  at  1500  °C  and  slow  cooling  to  room-temperature.  Whereas 
the  superstructure  reflections  characterizing  the  pyrochlore  long-range  ordering  of 
cations  and  anion  vacancies,  are  evident  in  the  first  ones  (e.g.  (1 1 1 )  or  (3  3 1 )  at  ~15° 
and  37°  (20)  respectively),  they  are  absent  in  the  latter  suggesting  increasing  struc¬ 
tural  disorder  as  Gd3+  atoms  are  replaced  by  the  smaller  Dy3+  ions.  Reported  XRD 
patters  of  F-  and  P-type  Gd2Zr207  as  well  as  that  of  our  pristine  Gd2Zr207  powders 
are  also  shown  as  a  reference.  Numbers  in  parenthesis  are  the  fluorite  and  pyrochlore 
Miller  indexes  of  each  reflection.  Left  and  right  inserts  show  an  enlargement  of  some 
areas  of  the  XRD  patterns. 


pyrochlore-type  zirconates  (Ln  =  Nd-Gd)  disorder  and  become  flu¬ 
orites  at  high  temperature  (e.g.  ~1550  °C  for  P-Gd2Zr207 ).  However, 
it  has  been  shown  that  even  in  zirconates  which  are  nominally  flu¬ 
orites  such  as  Y2Zr207  or  the  high  temperature  form  of  Gd2Zr207, 
random  mixing  of  cations  and  anion  vacancies  in  their  correspond¬ 
ing  substructures  is  never  complete  [14,15].  Therefore,  they  shall 
not  be  considered  as  “true”  fluorite-type  materials  but  better  rep¬ 
resented  by  ordered  pyrochlore-type  microdomains  in  a  disordered 
fluorite-type  matrix. 

As  shown  in  Fig.  1,  replacing  Gd  in  P-Gd2Zr207  by  homova¬ 
lent  Ln3+  cations  produced  either  anion  deficient  fluorites  or 
pyrochlores.  Thus,  Gd  substitution  by  smaller  lanthanides  such  as 
Dy,  Y  or  Er  (RGd  =  1.053  A  >RDy  =  1.027  A  >  RY  =  1.019  A  >  REr  =  1.004  A 
[16])  rendered  single  phase  materials  with  XRD  patterns  resem¬ 
bling  the  characteristic  of  anion  deficient  fluorites.  However, 
when  Gd  was  replaced  by  larger  cations  such  as  Sm,  Nd  or 
La  (RSm  =  1.079  A  <RNd  =  1.109  A<f?La  =  1.16  A),  the  XRD  patterns 
contained  all  the  superstructure  reflections  characteristic  of 
pyrochlores  (e.g.  the  (111)  and  (331)  reflections  at  ~15°  and 
38°  {20)  respectively).  Therefore,  decreasing  the  average  size  of 
cations  at  the  A-site  and  thus  the  value  of  the  RA/Rzr  ratio,  increases 
the  intrinsic  structural  disorder  of  P-Gd2Zr207  and  induces  a 
pyrochlore-to-fluorite  phase  transition.  In  contrast,  better  ordered 
pyrochlores  are  obtained  when  Gd3+  is  replaced  by  larger  lan¬ 
thanides  (the  Ra/Rzt  ratio  increases).  It  is  worth  mentioning  that 
although  cation  size  is  not  the  only  parameter  of  significance 
in  determining  the  degree  of  structural  ordering/disordering  in 
pyrochlores,  it  is  certainly  the  most  simple  and  convenient  for  the 
purposes  of  this  work. 

Peak  shifts  towards  higher  (Gd2_yDyyZr207)  and  lower 
(Gd2_ySmyZr207)  angles  observed  in  the  XRD  patterns  presented 
in  Fig.  1,  confirm  the  formation  of  the  title  solid  solutions  and  the 
expected  reduction  or  enlargement  of  lattice  parameters  when  Gd 
is  replaced  by  smaller  (Er,  Y  and  Dy)  or  larger  cations  (Sm,  Nd  and 
La)  respectively. 

As  for  the  electrical  properties,  the  Gd2_yLnyZr207  samples 
gave  rise  to  impedance  spectra  similar  to  those  typical  of  most 
solid  electrolyte  systems.  Fig.  2  shows  the  complex  impedance  plot 
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Fig.  2.  Complex  impedance  plot  obtained  at  300  and  350 °C  (inset)  for 
Gdi.7Ndo.3Zr207.  Solid  circles  represent  selected  frequencies  in  Hertzs. 

obtained  at  300  and  350  °C  (inset)  for  the  Gdi.7Nd0.3Zr2O7  com¬ 
position  where  two  main  features  are  evident,  i.e.  an  incomplete 
arc  at  low  frequencies  with  capacitance  values  typical  of  grain 
boundary  contributions  (~7  x  10-10  Fern-1  at  350 °C)  and  another 
one  at  high  frequencies  associated  with  the  bulk  contribution 
(~5  x  10_12Fcm_1  at  both  temperatures).  Bulk  dc  conductivities 
for  the  whole  series  were  calculated  from  these  experimental 
data  and  their  temperature  dependence  analyzed  by  using  an 
Arrhenius-type  law  of  the  form  crdc  =  (cr0/T)exp(-£dc//<r),  where 
£dc  is  the  activation  energy  for  the  conduction  process  and  cr0 
the  pre-exponential  factor.  Fig.  3  shows  such  representations  for 
selected  samples  where  the  solid  lines  are  least  squares  best  fits 
to  an  Arrhenius  law  confirming  that  ion  diffusion  in  these  mate¬ 
rials  is  thermally  activated.  The  dc  activation  energies  calculated 
from  the  slope  of  these  plots  were  found  to  vary  significantly 
with  Gd  substitution;  thus,  the  ^dc  value  of  1.13  eV  obtained 
for  undoped  P-Gd2Zr207,  increases  to  ~1.3eV  in  fluorite-type 
Gdi.2Er0.8Zr2O7  or  Gdi.2Dy0.8Zr2O7  and  decreases  to  1.0  eV  in 
pyrochlore-type  Gdi.6La0.4Zr2O7.  In  order  to  analyze  the  effect 
of  the  existing  structural  ordering/disordering  on  the  electrical 
properties  of  P-Gd2Zr207,  we  plotted  in  Fig.  4  the  calculated  dc 
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Fig.  3.  Arrhenius  plots  of  dc  conductivity  for  selected  samples:  Gdi.2Ln0.8Zr2O7 
(Ln  =  Er,  Dy,  Sm,  Nd,  La),  Gd2Zr207  and  Gdi.6Yo.4Zr207.  Solid  lines  are  the  least  squares 
best  linear  fits  to  the  experimental  data. 


RA/RGC 

Fig.  4.  dc  activation  energies  for  oxygen  ion  diffusion  in  Gd2-yLnyZr207  plotted  as  a 
function  of  a  normalized  ionic  radius  of  the  A  cation  where  Ra/Rga  =  1  corresponds 
to  the  undoped  material.  The  dashed  line  is  not  a  linear  fit  and  is  only  shown  to 
emphasize  a  trend. 


activation  energies  for  the  whole  series  as  a  function  of  a  nor¬ 
malized  average  A  cation  ionic  radius.  To  calculate  the  average  RA 
size  in  each  sample,  we  used  the  individual  ionic  radii  given  by 
Shannon  [16]  as  well  as  the  fraction  of  each  cation  present:  e.g., 
the  average  RA  size  in  Gdi.7Nd0.3Zr2O7  would  be  calculated  as  fol¬ 
lows:  Ra  =  [  1.7i?cd  +  0.3KNd]/2.  To  facilitate  the  estimation,  we  have 
assumed  no  mixing  between  cations  at  the  A  and  B-sites  that  is, 
Gd3+  and  Ln3+  ions  are  constraint  to  the  A-site  whereas  Zr4+  is  exclu¬ 
sively  located  at  the  B  position.  Finally,  to  ease  the  discussion  we 
have  also  normalized  the  RA  value  in  each  sample  relative  to  that 
of  Gd3+  and  plotted  the  experimental  data  versus  the  RAIRGd  ratio. 

As  Fig.  4  shows,  the  Edc  needed  for  oxygen  migration  increases 
significantly  with  decreasing  RAlRGd  (increasing  structural  disor¬ 
dering);  in  fact,  the  Edc  follows  a  general  trend  which  persists 
even  across  the  theoretical  fluorite/pyrochlore  phase  boundary  in 
the  system,  and  shows  a  sharp  and  almost  linear  decrease  with 
increasing  RAIRGd.  As  cell  volume  in  fluorite  and  pyrochlore-type 
Ln2Zr207  zirconates  is  to  a  good  approximation,  linearly  depen¬ 
dent  on  the  ionic  radius  of  the  A-type  cation  [17],  Gd  substitutions 
in  P-Gd2Zr207  by  larger  cations  will  increase  cell  volume  and  thus, 
facilitate  ion  hopping.  Therefore,  a  lower  energy  barrier  for  mobile 
oxygen  ions  would  be  expected  which  is  consistent  with  the  exper¬ 
imental  data  of  the  title  series  and  with  those  published  previously 
in  similar  pyrochlore  systems  [18]. 

Fig.  5  shows  the  evolution  of  the  bulk  conductivity  in  the 
Gd2_yLnyZr207  system  versus  the  normalized  RAIRGd  ratio  at  two 
measuring  temperatures,  500  and  800  °C.  As  this  figure  shows,  con¬ 
ductivity  in  the  series  does  not  depend  only  on  the  dopants  size 
and  their  influence  in  cell  volume  but  also  in  their  effect  on  oxygen 
vacancies  ordering.  Thus,  conductivity  at  any  given  temperature  is 
higher  for  the  pyrochlore-type  materials  than  for  those  showing 
the  anion  deficient  fluorite  structure.  Within  each  structural  type 
however,  conductivity  values  at  high  temperature  remain  almost 
constant  and  independent  of  the  dopant  cation  and  composition. 
As  Fig.  6a  and  b  shows,  these  changes  in  conductivity  cannot  be 
explained  by  differences  in  microstructure  since  all  the  samples 
analyzed  showed  after  sintering,  similar  final  porosity  and  grain 
size. 

Finally,  we  have  not  observed  a  peak  in  conductivity  close  to 
the  pyrochlore/fluorite  phase  boundary  as  suggested  by  some  other 
works  carried  out  in  similar  Ln2Zr207  systems  [17]  although  con- 
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Fig.  5.  The  effect  of  homovalent  Gd  substitution  in  the  bulk  conductivity  of  the 
pyrochlore-type  Gd2Zr207,  plotted  as  a  function  of  a  normalized  averaged  A  cation 
size. 

ductivity  in  the  first  is  certainly  higher  than  in  the  latter.  Many 
papers  have  been  published  dealing  with  the  effect  of  dopant  size 
on  the  conductivity  of  different  ionic  conductors  and  it  is  gener¬ 
ally  accepted  that  those  inducing  the  lowest  strain  in  the  crystal 
lattice  of  the  host  material,  exhibit  the  highest  conductivity.  Thus, 


Fig.  6.  Scanning  electron  micrographs  of  Gd2Zr207  (a)  and  Gdi.6Lao.4Zr207  (b),  after 
pressing  and  sintering  as  described  in  Section  2. 


in  the  (Gdi_xAx)2Ti207  system  (with  A  =  Mg2+,  Ca2+,  Sr2+,  K+),  the 
highest  conductivity  values  were  obtained  for  Ca2+,  which  exhibits 
the  closest  size  match  to  Gd3+  [19].  Similar  studies  carried  out  in 
fluorite-type  Lny  Cei  _y  02_«$  have  also  shown  a  minimum  in  the  acti¬ 
vation  energy  and  a  maximum  in  conductivity  for  Gd3+,  with  the 
closest  size  match  to  Ce4+  [20].  In  our  case,  the  pristine  material  P- 
Gd2Zr207  is  not  an  insulator  as  Gd2Ti207  neither  doping  is  used  to 
create  vacancies  as  in  Ce02  but  to  manipulate  the  degree  of  struc¬ 
tural  ordering/disordering.  Our  results  show  that  P-Gd2Zr207  is 
highly  tolerant  to  chemical  substitutions  at  the  A-site  with  the  best 
dopants  being  those  which  increase  the  size  mismatch  between 
cations  at  the  A  and  B-sites  that  is,  lanthanide  ions  larger  than 
Gd3+.  Decreasing  number  of  mobile  charge  carriers  with  increas¬ 
ing  structural  ordering  is  apparently  compensated  by  decreasing 
activation  energies  for  migration  and  thus,  conductivity  remains 
almost  constant  (or  even  increases  slightly)  as  observed  in  the 
Gd2(Tii_yZry)207  system  when  Zr4+  is  replaced  by  Ti4+  [21,22]. 

4.  Conclusion 

We  have  shown  that  pyrochlore-type  Gd2Zr207  is  highly  toler¬ 
ant  to  a  large  number  of  A-site  lanthanide  substitutions.  Whereas 
replacing  Gd  by  smaller  lanthanides  yields  anion  deficient  fluorites, 
larger  cations  produce  increasingly  ordered  pyrochlores.  Interest¬ 
ingly,  the  Gd2_yLnyZr207  solid  solutions  (Ln  =  Sm,  Nd  and  La)  are 
good  ionic  conductors  with  smaller  dc  activation  energy  for  oxygen 
migration  than  that  observed  for  the  pristine  P-Gd2Zr207. 
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